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Abstract 

Data on J/\I/ — > 7(?77r + 7r _ ) and 7(A =t Ag7T =F ) from 58M J/4' hadronic interactions in 
the BES II detector are analysed. They throw new light on j p = o- mesons throughout the 
mass range up to 2 GeV, notably r/(1440). For the first time, this state is fitted including 
the full s-dependence of decays to ao(980)7r, rja, KK*, kK and /o(980)r/, including the 
dispersive term associated with this s-dependence. Two types of fit are reported. The 
first uses a single ?y(1440). The second uses separate 7/(1405) and 7/(1475). The strong 
s-dependence of P-wave production of K*K shifts the fitted mass of 7/(1475) down to 1440 
MeV. As a result, the two types of fit give almost identical results and the conclusion from 
s * ! these data is that there is no significant evidence for two separate resonances. In KKir, 

there are definite though small contributions from /i(1420) and /i(1285). The optimum 
fit requires either a small additional /i(1510) or alternatively a small perturbation to the 
high mass tail of /i(1420). At higher masses, there is a broad J p = 0~ peak in 7/7T7r at 
1850 MeV. There is also a weak but definite signal consistent with 772 (1870) — > 7/7T7T. 

i> : 

43 , . PACS: 13.25.-k, 14.40.Cs 

^ ! 1 Introduction 

The analysis reported here was done in the years 2002-3 in collaboration with two senior mem- 
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bers of the BES collaboration. The work was funded by the Royal Society and Queen Mary 
College via an official agreement with the Chinese Academy of Sciences and the BES collabo- 
ration. This agreement guaranteed access to BES 2 data for the purpose of publications. Since 
then, the BES management has refused to appoint internal referees, hence blocking a collabora- 
tion publication and comments from the collaboration. The data illuminate the spectroscopy of 
J p = 0~ mesons and should be in the public domain. There is a responsibility to publish work 
which has been supported from public funds with approximately £60, 000 . The wording of the 
article is close to a document submitted late in 2008 to the BES management, except for some 
clarification of wording. 

There have been many earlier studies of r/(1440) in J/4 7 radiative decays, pp annihilation 
and 7rN interactions. Extensive references are to be found in the listings of the Particle Data 
Group (PDG) [1]. Their review concludes that two nearby 0~ + resonances are required in the 
mass range 1400-1500 MeV. The lower one is at 1405 MeV, decaying to rj-n-ii and kK, where k 
stands for the Ktc S-wave. The upper state is at 1475 MeV, decaying dominantly to A"* (890) A". 

However, almost all of these analyses are based on fitting the strong A* (890) A channel 
with a Breit-Wigner resonance of constant width. The threshold for this channel is at 1390 
MeV, although smeared out to some extent by the 50 MeV full- width of the A*. The A* A 
channel with J p = 0~ is produced in a P wave. As a result, the intensity of the A*A channel 
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increases with kaon momentum k as k 3 . This rapid s-dependence enhances the upper side of 
the K*K peak strongly. It also gives rise to a strong dispersive correction in the Breit-Wigner 
denominator, as described in subsection 4.3. Only one previous analysis includes the first feature 
and none contains the second. They are therefore prejudiced towards two separate resonances. 

As well as the possible 7/(1405) and 77(1475), there is evidence from four experiments on 
up — > (i]7ttt)p for 77(1295) [2], [3], [4] and [5]. In these data, the possibile exchange of both natural 
and unnatural parity mesons requires a rank 2 analysis, with partial interference between the 
two alternatives, hence incomplete interferences between /i(1285) and 7/(1295). The possibility 
of no interference between /i(1285) and 77(1295) leaves some room for uncertainty in conclusions. 
There is no definite evidence for 77(1295) in J/^J decays and only weak evidence mpp annihilation 
[6]- 

There are therefore three controversial candidates for two radial excitations of 77 and 77'. 
Some people have argued that 7/(1295) has been confused with /i(1285) [7]. Others argue from 
the strong production of 7/(1405) and 7/(1475) in J/\l/ radiative decays that one of them is the 
0~ glueball [8]. However, Lattice Gauge calculations predict the glueball above 2 GeV [9]. 

One must remember too that the 77 and 7/ are well established to have compositions 

77 = nn cos — ss sin (1) 
77' = 7777 sin + SS COS 0, (2) 

where 7777 = (uu + dd)/V2 and = 41.5° if there is no glueball component [10]. In J/ty radiative 
decays, the 77' is produced with an intensity a factor 5 larger than the 77. Similar mixing for 
radial excitations could explain surprisingly strong production of 7/(1405) and/or 7/(1475) in J/\I/ 
decays and pp annihilation. 

There are two likely scenarios. The first is that the 7/(1295) does not exist and 7/(1405) may 
be assigned dominantly as 7777 and 7/(1475) as ss. The second scenario is that 7/(1295) does exist 
and 77(1405) and 7/(1475) are really one resonance. These two alternatives are fitted to data 
presented here. The s-dependence of the Breit-Wigner amplitudes is treated fully for the first 
time. 

The higher mass range for i]tttt is also of considerable interest. BES II data provide evidence 
for a 7 standard deviation signal in rfTrir for a narrow X(1835) with T = 68 MeV [11]. It is 
important to know if it appears in i]7ttt and if so with what branching fraction. 

Section 2 presents technical details of the selection of data for the channels 7(7/7r+7r~), 
r )(nK + K^) and ^(^K^tt^) and discusses backgrounds. The r )(nK + K~) channel has not been 
examined before. It is useful in establishing the magnitude of possible /o(980) contributions to 
rjKK. Sections 3 and 4 present features of the data and the formulae needed to treat the s- 
dependence of Breit-Wigner amplitudes. Section 5 gives details of the fit to the mass range up to 
1600 MeV. This is the centre-piece of the paper. Section 6 considers the high mass range above 
1600 MeV and estimates the magnitude possible for V(1835) — > r/mr; the low value observed is 
inconsistent with the X(1835) being dominantly 7777. Section 7 returns to the interpretation of 
7/(1440) and /1 (1510). Section 8 summarises conclusions. 
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2 Data selection 



The data fitted here were taken with the BES II detector. Full details of the detector and its 
upgrade are reported by Bai et al. [12], [13]. It has cylindrical symmetry around the intersecting 
e + e~ beams. Its essential features are (i) a Main Drift Chamber (MDC) for the measurement of 
charged particles, (ii) time of flight detectors with a o of 180 ps, and (iii) a 12 radiation length 
Barrel Shower Counter comprised of gas proportional tubes interleaved with lead sheets. The 
MDC measures dE/dx. Together with the time of flight detectors, it separates ir and K up 
to ~ 700 MeV/c. Outside the 7 detectors is a magnet providing a field of 0.4T. This magnet 
is instrumented with muon detectors, but for present work they serve simply to reject 
pairs. The Main Drift Chamber provides full coverage of charged particles for lab angles 9 with 
|cos#| < 0.84. 

2.1 Selection of ^(r]7T + 7T^) and ~f{r}K + K~) 

These events have been selected by demanding exactly three photons and two charged tracks 
with balancing charges. Both charged particles must be positively identified by time-of-flight 
and/or dE/dx as pions (for r]n + iT~) or as kaons (for r/K + K~). The vertex must lie within a 
cylinder of 1.5 cm radius, 30 cm long, centred on the intersecting e + e~ beams. Two photons 
must have an invariant mass within 70 MeV of the 77. Events are rejected if any two photons 
produce a mass within 40 MeV of the 7r°. Photons may also originate from interaction of charged 
tracks in the detector. Any photons lying within a 6° cone around a charged track are discarded. 

Candidate events have been subjected to 4-constraint kinematic fits to 377r + 7r _ and 5C 
fits to 7?77r + 7r~ (or correspondingly 3 r yK + K' and r yr]K + K^); they are accepted only if the x 2 
probability of each fit is > 5%. Any events where two photons are consistent with the final state 
2 / yr]7r + 7T~ or 2 / ~fi]K + K~ are rejected if x 2 < 20. They are also rejected if the measured particles 
are consistent with the loss of one photon and a fit to 477r + 7r~ or 4 / yK + K~ with x 2 < 20. 

The main background arises from the loss of one photon from the final state rj-n + -n'-n° , which 
is dominated by a 2 (1320)p. This background has been studied using the standard SIMBES 
Monte Carlo simulation of the detector. After the loss of one soft photon, it leads to a strong 
peak at high i]7i7i masses 2300-2800 MeV. This background may be reconstructed from data by 
assuming one missing photon in the kinematic fit. For events fitted this way, Fig. 1(a) displays 
the ir ^ mass distribution. There is a clear peak due to p(770). If M(7r°7r ± ) lies in the mass 
range 670-870 MeV, one finds the 7771- T mass distribution of Fig. 1(b). It shows clear a (980) 
and a2(1320) peaks. In order to suppress these backgrounds further, the missing momentum is 
calculated from the measured i]7t + tt~ of each candidate event. Assuming the missing particle 
to be a 7T°, any event with M(7r°7r ± ) in the mass range 670-870 MeV is rejected. However, a 
minor problem is that events lying in the charge zero peak of the p cannot be eliminated without 
rejecting required 7770- events, where a is a shortand for the nn S-wave. 

Fig. 1(c) shows the mass distribution for the 77 pair closest to the 77 for all combinations 
fitting 377r + 7r~. Fig. 1(e) shows the corresponding plot for candidate r]K + K~ events. Fig. 1(d) 
shows the rjnn mass distribution from the final data sample. Peaks are visible at 1285 and 1415 
MeV, followed by a sharp rise at 1490 MeV to a shoulder at ~ 1550 MeV. There is a further 
peak at ~ 1850 MeV, then a rapid fall to 2 GeV. 

Above 2 GeV, there is a large contamination from the background explained above. The 
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Figure 1: Mass distributions for (a) and (b) r]7i T for events fitted to 7r°^7r + 7r _ after 

selecting events with M(n 7r ± ) in the range 0.67 to 0.87 GeV; (c) mass distributions for the 
77 pair closest to the rj for 377r + 7r~ and (e) "i^K+K'; (d) M(r]-K + ir~) after a kinematic fit 
to / yi]tt + 7t~; (f) M(i]K + K~) after a kinematic fit to / yi]K + K~; the dashed curve shows the 
predicted mass distribution for / (980) — > K + K~ from ?y(1440). 



Monte Carlo simulation of the a 2 p background is fitted to events in the mass region 2300- 
2800 MeV, hence establishing its magnitude. It agrees with the PDG branching fraction of 
J/\l/ — ► a 2 (1320)p within its quoted error of 35%. 

The physics analysis will be restricted to the mass range below 2 GeV. Since a 2 (1320)p is the 
dominant background by a factor 10, the background at low r]n + Ti~ masses is taken from the 
Monte Carlo simulation of this channel. Integrated up to an i]7iiT mass of 2 GeV, this background 
is 18.2%. Its distribution with r/mr mass is illustrated below in Fig. 7(d). In the mass range 
below 1600 MeV it is small and featureless. Although there is a background of misidentified 
events, direct feedthrough from a 2 p produces at most 33 events in p — > 7r + 7r~, mostly above 
1600 MeV. This small a 2 (1320) background surviving in data will be pointed out later. 

At the end of the data selection, there are 18,751 events up to an rjinr mass of 2000 MeV 
of which an estimated 3350 are background. The average efficiency for data reconstruction is 
13.5%. The branching fraction is 

BR[J/m -> 7(^7T7r)] = (7.5 ± 0.8(syst)) x 10~ 3 (3) 
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up to a mass of 2 GeV, after correction for all charge combinations and background. The error 
is mostly systematic. It covers the estimated background and the uncertainty in partial wave 
amplitudes due to loss of events through the 35° holes in the detector around the beam entrance 
and exit. The PDG quotes an average branching fraction (6.1 ± 1.0) x 10~ 3 . However, this refers 
to a Crystal Ball measurement [14] claiming a peak in r/nn at 1700 MeV, substantially different 
to what is observed here. 

2.2 Data on J/V -> ^{t]K + R-) 

For this final state, the main background arises from 7r°i]K + K~ after the loss of one photon. 
The contamination is estimated to be ~ 27% from a study of side-bins either side of the r) 
peak, together with a Monte Carlo estimate of the shape of the background. The reconstruction 
efficiency is 7.6% and varies little with r/KK mass. Fig. 1(f) shows the mass distribution of 
events fitting as ^r]K + K~. The number of events is too small to allow a meaningful physics 
analysis. However, the data are useful in establishing an upper limit on the number of events due 
to the final state /o(980)t/, /o — » K + K~ . The essential point is that the number of signal events 
is a factor 100 smaller than for 'yrjii^Ti^. After background subtraction, there are 22 events 
up to an r/KK mass of 1.8 GeV and 125 up to 2 GeV. The latter corresponds to a branching 
fraction 

BR{J/^ -> ~fr]KK} = (1.4 ± 0.4) x 10" 4 (4) 
up to 2 GeV. Errors are again systematic, from the acceptance. 

2.3 Selection of 7 (A' ± A'§7r T ) events 

These events have been selected initially with at least one photon and with 4 charged tracks 
having balancing charges. One and only one charged kaon must be identified by time-of-flight 
or dE/dx; at least two charged tracks must be positively identified as pions. 

Events are then subjected to a AC kinematic fit, requiring a x 2 probability > 5%. The mass 
of one 7r + 7r~ pair must lie within 25 MeV of the K s . This pair must be consistent with a vertex 
within 6 cm of the second vertex defined by the other two charged particles. Fig. 2 shows the 
mass distribution for the reconstructed Kg before the kinematic fit. There is a clean peak with 
only a 4% background, which can be evaluated accurately. Events are rejected if the spectator 
pion and either of the other two particles (fitted as pions) have an invariant mass within 25 
MeV of the K° s . 

Events are rejected if there is a better fit to final states K ± K s tt t or K ± K s -k t tt° without an 
additional photon. In order to reject the second process more strongly, events are discarded if 
the missing mass attributed to neutral particles is within 50 MeV of the 7r° mass. 

The kinematic fit is made assuming that the K s decays at the primary vertex; the effect 
of this assumption on the selection of events is simulated by a Monte Carlo study. There is 
a small loss of events, but no significant bias to kinematics. Alternative selections have been 
made fitting two separate vertices. These lead to a larger loss of events because of the difficulties 
of reconciling two vertices with the kinematics; within statistics, there is no change to physics 
conclusions. 

At the end of the data selection, there are 5638 events including an estimated background 
of 1068. The background is substantially smaller over the reduced mass range up to 1600 MeV. 
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Figure 2: The mass distribution for the reconstructed K s . 



The average efficiency for data reconstruction is 10.8%. The resulting branching ratio, corrected 
for charge combinations, is 



BR[J/m -> jKKir] = [2.2 ± O.S(syst)} x 10" 



(5) 



up to a mass of 1800 MeV. The error is almost entirely systematic. A 10.3% error arises from 
extrapolating the partial wave analysis over the holes of the detector, 7.5% uncertainty in the 
effects of cuts, and 4.5% in the background; these errors are added in quadrature. The PDG 
quotes an average of (2.8 ± 0.6) x 10~ 3 , though this is derived from four values spanning the 
range (1.66-4.3) x 10~ 3 . 
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Figure 3: The K ± K S K T mass spectrum. 



The Monte Carlo used here to simulate the BES detector was the standard SIMBES package 
released in February 2002, in its version as of March 2003. During discussions with the BES 
management, they have claimed that alterations to the calibration of the time-of- flight system 
after March 2003 could affect the selection of events. This point is discussed in detail in the 
accompanying paper on 7(7r + 7r~7r + 7r~) data [15]. In that analysis, no systematic discrepancy 
is observed between (dE/dx) and time-of-ffight, and the same is true for present data; no 
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contamination of 77r + 7r~7r + 7r° events by r yK + K~n + n~ is observed above the 1% level. In the 
present data, there is again no observed confusion between 7r + 7r~ and K + K~ . For jr]7i + n^ 
data, contamination from ^r]K + K~ is obviously < 1% because of the small observed branching 
fraction for that channel. There is no evidence for contamination from ir°i]K + K~ , which is 
kinematically well separated from / yi]TT + TT^; such contamination is estimated to be < 1%. For 
/ j(K ± K s it t ) data, the selection of the K s limits possible background to 4%. The observed 
background under the K s is consistent with that originating from tt° K ± K s n T ; it rises with 
KKir mass as expected for that source. 

3 Features of the data 

Fig. 3 shows the KKn mass distribution in J/ty — > r y(K ± K s 7i T ). A small f\ peak is visible 
at 1285 MeV. Then there is a strong but asymmetric peak at ~ 1450 MeV on a slowly rising 
background. There may be a small structure causing a drop at 1530 MeV, though this could be 
a large statistical fluctuation. A similar shoulder is visible at 1510 MeV in Mark III data of Bai 
et al. [16]. Several alternative prescriptions for fitting this shoulder will be described. 

The raw mass distribution for r\ix+ix~ is shown in Fig. 1(d) and fits to it will be shown in 
detail below in Figs. 7 and 11. There is a narrow peak at 1415 ± 10 MeV with a full width of 
60 MeV. There is also a sharp rise at 1500 MeV. It is vital to realise from the outset that these 
structures are superimposed on a large and broad 0~ signal extending over the whole mass range 
from 1200 to 2000 MeV and interfering coherently, see Fig. 7(b) below. Without the inclusion 
of this broad component, fits fail completely to describe the rjnn mass projection. In particular, 
the broad component plays a major role in fitting the shoulder at 1550 MeV. 

The 747r data contain a similar broad peak at 1600 MeV in pp with J p = 0~ [15]. It is 
natural that some pp pairs will de-excite to r/cr or a (980)7r with L = 1 in the decay. It is also 
possible that the pp peak at 1600 MeV is resonant. 

Dalitz plots for r]iT + iT~ are shown for three mass ranges in Fig. 4. The upper 3 panels are 
for data and the second three from the fit to 7^(1440) reported below. The lower set does not 
contain the statistical fluctuations of the upper set and the plotting procedure makes it appear 
that there are some differences between data and fit. However, a check on all bins shows no 
evidence for systematic discrepancies. The main feature in all three cases is a peak due to two 
crossing a (980) bands. 

Dalitz plots for KK-rr are shown in Fig 5 for three mass ranges. Here, the main feature is due 
to two crossing i^*(890) bands. There is a strong enhancement near the intersection of these 
bands in Figs. 5(a) and (b). It is much reduced in Fig. 5(c). This is due to the zero in the Kir 
amplitude near the middle of the K* band, arising from the angular dependence of the partial 
wave amplitude in the decay. 

4 Partial wave analysis 
4.1 Preliminaries 

The aim of the analysis is to fit the data in terms of conventional resonances, plus the broad 
interfering J p = 0~ component observed both here and in the 7(7r + 7r~7r + 7r~) data. It would be 
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Figure 4: Dalitz plots for r/mi: (a)-(c) show data and (d)-(f) the fit. Mass ranges for rjirir are 
1350-1480 MeV in (a) and (d), 1480-1620 MeV in (b) and (e), 1620-1950 MeV in (c) and (f). 



good if the spin composition of rjmr and KKix could be determined in narrow slices of mass, 
and this has been tried. The main conclusion from that attempt is that 0~ is strongly dominant 
for both sets of data. In KKix data, contributions from J p = 1 + are identified clearly in the 
narrow mass ranges of /x(1285) and / 1 (1420). There is also definite evidence for a small 1 + 
contribution continuing above /i(1420), up to ~ 1550 MeV. However, it is sufficiently small 
that its line-shape cannot be determined with any precision from slice fits. The contribution 
from J p = 2~ is small. 

There are however two difficulties with trying to analyse the data in slices. The first is that 
the only available information is a differential cross-section. This does not allow independent 
determinations of amplitude and phase. A related difficulty is that there are complex interfer- 
ences between ao(980)7r, rja and f (980)i] in -qmr data and likewise in KKix data between K*K, 
kK and ao(980)7r, ao(980) — > KK. In the absence of any strong interferometer in the data with 
known quantum numbers, the slice analysis cannot sort out these interferences. This leaves no 
choice but to relate magnitudes and phases through Breit-Wigner amplitudes fitted directly to 
the data. 

Partial waves which are tried are J p = 0~, (orbital angular momentum V — 1 in the 
production process J/V -> 7 + X), 1+ (V = 0), 2~ (V = 1) and 2+ (V = 0). Of these, J p = 2+ 
is negligible. Decay channels which are considered are ao(980)7r, K*(890)K, kK, a2(1320)7r, 
/ 2 (1270>7, / (980>7 and W . 

There is an immediate question how a and k should be parametrised. In nn and Kn elastic 
scattering, there are Adler zeros which suppress the amplitudes strongly near threshold. This is 
well understood in terms of Chiral Symmetry Breaking, which creates the Adler zeros. However, 
in many sets of data, it is established that a and k, appear as low mass peaks. The origin of this 
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Figure 5: Dalitz plots for K^Kgir^: (a)-(c) show data and (d)-(f) the fit. Mass ranges for 
KKn are 1380-1440 MeV in (a) and (d), 1440-1500 MeV in (b) and (e), 1500-1560 MeV in (c) 
and (f). 



effect is straightforward. If the elastic amplitude is written as N(s)/D(s), the denominator must 
be universal. However, the numerator can be quite different in elastic scattering and production 
reactions. BES II data for J/ty — > luit + it~ [17] may be fitted accurately with N(s) a constant. 
This is also true for BES data on -> KKtttt, where there is again a low mass peak in the 
Kit S-wave due to the k pole [18]. The present work uses the k amplitude fitted in Ref. [19] to 
LASS data for elastic scattering [20], BES II, and E791 data. 

A problem is that it is not yet fully understood why or how the Adler zero disappears 
from the numerator of the amplitude in some production processes. In principle, there could be 
contributions with a numerator N(s) = A + Bs where A and B are complex constants. However, 
this unfortunately leads to large destructive interferences between the two terms. This is a well 
known symptom of over-fitting the data, so final fits must choose between a and k poles or 
alternatively amplitudes the same as for elastic scattering. Empirically, the present ^(KKit) 
data give a considerably better fit (by 46 in log likelihood) with the k pole amplitude, i.e. a 
constant numerator N(s). [Log likelihood is defined here so that it changes by 0.5 for a one 
standard deviation change in the fit]. The 7(?77r + 7r~) data show less discrimination, but still 
prefer a fit with the o pole, i.e. no Adler zero in N(s). This fit is better than one using the 7T7t 
elastic scattering amplitude by 16 in log likelihood. Fits to present data are almost identical 
using the formula for the a of Zou and Bugg [21] or a recent formula fitted to the a pole [22]. 
The first formula is simpler and is used here. Both a and k are fitted using a constant for N(s). 

The fit to data is made using 200K Monte Carlo events for r )(rj'n + 'K~) and 52K for ^(^Kg^) 
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to simulate detector acceptance. These Monte Carlo events are generated with the standard 
SIMBES package. All figures show experimental results uncorrected for acceptance. These cor- 
rections are included automatically into the log likelihood fit. The acceptance varies only slowly 
across figures. 

A technical detail is that partial wave amplitudes are written in terms of relativisic tensors; 
explicit formulae are given in [23]. In J/\l/ radiative decays, an important point is that the 
dominant production process is expected to be via intermediate •ycc. Because of the high mass 
of these quarks and the interaction of the photon at a point, the production process is expected 
to be almost pointlike. Data on J/$ ^ ^An in the accompanying paper verify this [15]. They 
contain a dominant J p = CT signal in pp with a dependence on An mass close to the expected 
P 3 for 0~ production with V — 1; here P 7 is the photon lab momentum. The consequence 
for present data is that production should proceed only through the lowest available orbital 
angular momentum V . There are well identified /i(1285) and /i(1420) signals in ^(KKtx) final 
states. In both cases, any amplitudes with V = 2 are negligible, producing improvements in log 
likelihood < 4. For production of 0~ states, there is likewise no evidence for L' = 3. With this 
simplification, there are distinctive differences between final states with J p = CT, 1 + , 2~ and 
2+ 



4.2 Details of partial wave fits 

Fits to the 1450 MeV KKir peak immediately reveal the need for both J p = 1 + and 0~ 
amplitudes. The /i(1420) amplitude is conspicuous from its cos Od dependence, where 6d is the 
decay angle in the rest frame of the K* to K and n. A free fit gives M = 1429 ± 3 MeV, 
compared with the value 1426.4 ±0.9 MeV of the PDG. If the /i(1420) is removed from the final 
fit, log likelihood is worse by > 129 for two less fitting parameters. Statistically this is nearly a 
14 standard deviation effect. 

There is also a small but definite /i(1285) contribution, decaying to kK. If the 7/(1295) is 
included in the fit, its contribution is small and will be discussed in Section 7. 

Earlier analyses have freqently fitted the a (980) freely to KKir data. The a (980) contri- 
bution lies along the upper right-hand edge of Dalitz plots of Fig. 5. For low KKn masses, it 
interferes strongly with the two K*(890)K components, see Figs. 5(a) and (b). It can easily 
be confused with the intersection between the two if* (890) bands. The a (980)7r contribution 
in r/7T7T data does not suffer from this problem. Fits have therefore been made constraining the 
decay branching ratio of a (980) between r/ir and KK to a value consistent with the Flatte form 
for ao (980) [24]. Integrated over the ao(980) peak up to 1.3 GeV, the branching fraction to KK 
is (24 ±2)% of that in rjn. This constraint results in small but well defined contributions of 
a (980) to KKn data. Production phases are fitted freely. 



4.3 Parametrisation of amplitudes 

It is essential to include into the Breit-Wigner amplitude the strong s-dependence of the phase 
space for P-wave production of K*K. This phase spaces rises as k 3 , where k is the kaon 
momentum in the K*K rest frame. As a reminder, the intensity for a two-channel resonance is 

I(s) = MUM (6) 



10 



where T tot is the total width. The factor T 2 (s) in the numerator arises from the phase space of 
the final K*K state. The same s- dependence mustbe included into the denominator, otherwise 
an inconsistency with unitarity arises between numerator and denominator. The correct form 
for the Breit-Wigner amplitude for decay to channel j (for example K*K) is 

f = FjBjpjjs) 

h M 2_ s _ m(s) - lYjl gtp t {s)F?Br 1 ) 

where <?j is the coupling constant for each channel % and p(s) is its phase space. A form factor 
F is included, and (where required) a centrifugal barrier factor B. Details will be given shortly. 
The evaluation of phase space is done numerically, using Eq. (40) of Ref. [25] which will be 
repeated here for completeness. As an example, the 3-body phase space for K*K is given by 

Pk*k{s) = / „ ^^\T K * R {si)\ , (8) 



4m 2 K 7T y/SSi 

where T is the Breit-Wigner amplitude for the K*. Here s refers to the resonance X and k 
to the momentum of the K* or K in the rest frame of X; s\ and k\ refer to the K* and the 
momenta of decay K and n in its rest frame. For cases where there is angular momentum in 
the decay, the centrifugal barrier factor B appears in the numerator. 




1 1.2 1.4 1.6 1.8 



Mass (GeV) 

Figure 6: The dispersive term m(s) for ^(1440). 



The term m(s) is a dispersive correction to the real part of the amplitude, required to make 
the amplitude fully analytic: 

m(s) M2 - s f MT tot (s')ds> 

m[S) ~ 7T J (S> - S)(M* - S'Y W 

It includes a subtraction on resonance making the integral strongly convergent. The resulting 
m(s) for 77(1440) is illustrated in Fig. 6. In the 1440 MeV mass region, the term m(s) actually 
varies more strongly with s than the term (M 2 — s) in the Breit-Wigner denominator, leading 
to major deviations of phases from those of a Breit-Wigner resonance of constant width. This 
is one of the essential new features of the present analysis. 

In fitting rjim data, it is unavoidable to introduce a broad J p = 0~ signal. It needs to 
accomodate peaks observed in J/ty — > pp, luuj, K*K* and 00 near their thresholds. A natural 
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explanation is that it arises from J/\I/ — ■> 7GG, where G are gluons, followed by coupling of 
gluons to vector mesons by colour neutralisation. An essential feature of the 7(7^71-) data is 
the requirement for a peak at ~ 1600 MeV, naturally associated with \pp\ L= i — > [r]a] L= i and 
[a (980)7r] i= i. The strategy adopted here is to examine how far one can get by fitting this broad 
component by a single broad resonance X in the production process; exactly the same formula as 
(7) is used to fit 7(471") data in the accompanying paper [15]. The mass M optimises empirically 
at (2.04 ± 0.10) GeV. This is far removed from 77(1440) so the only significant phase variation 
arises from (important) dispersive effects associated with the opening of each threshold. This 
broad background will be described as 'the broad - '. 

Although Eq. (7) may look complicated, it is straightforward to programme and gives an 
excellent fit to present data and simultaneously to data on J/$ ^ 747T. Each of the channels r/cr, 
ao(980)7r, /o(980)7r, K*K, kK, pp, lulu, K*K* and 00 requires just one parameter g 2 constrained 
(iteratively) to reproduce observed branching fractions of these channels. This constrains the 
s-dependence of the broad 0~ amplitude closely. The first five channels are essential to fit rj-Kn 
and KKtt data. The inclusion of the K*K* threshold examines whether or not the drop in the 
i]7i7T data above 1850 MeV can be explained by the K*K* threshold. The channels lulu and 00 
have little effect but are included for completeness. 

Further details concern the form factor F and centrifugal barriers B in Eq. (7). The form 
factor for J/^> — > X, X — > j is the product of a weak form factor F p for the first step and a 
factor F d for decays of X through channel j. Both are taken as Gaussians: 

F = exp(—a p K 2 ) exp(— a^A; 2 ), (10) 

where K is the production momentum of X in the J/\l/ rest frame and k is the decay momentum 
in the rest frame of X (for example, the momentum of 77 in the channel 770"). Values of a p and a<j 
are taken from the work on J/ty — > 747r: a p = 0.06 GeV -2 and a d = 2.25 GeV -2 , corresponding 
to RMS radii R p = 0.12 fm, R d = 0.73 fm. The production form factor is very weak and has 
little effect here, but is included for consistency with 7(477) data. 

The factor B is likewise the product of the L' = 1 centrifugal barrier for J/^> — > 7X and 
(where necessary) a centrifugal barrier for the decay, for example an L = 1 barrier for the decay 
of 0~ to K*K. As a reminder, for L = 1, 

B = kR/Vl + k 2 R 2 - (11) 

the same value of R is used in the centrifugal barrier as in the form factor, though this may be 
an approximation. 

5 Fits to data up to 1600 MeV 

Two types of fit will be discussed. The first (Fit A in tables) tries to fit both 777T7r and KKir data 
with a single 7/(1440) resonance. The second (fit B) uses two separate resonances like 7/(1405) 
and 7/(1475), though their masses and widths need re-optimising. 

Fig. 7 shows details of the fit to ^(r/Trn) data with a single resonance which optimises at 
a mass of 1439 ± 5 MeV. Panel (a) shows the overall fit to the mass distribution; the lower 
histogram shows the coherent sum of 7/(1440) and the broad Panel (b) shows the full 0" 
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Figure 7: The fit to ?77r7r data using r/( 1440) . Points with errors show data. In (a) the upper 
histogram shows the full fit; the lower one shows the coherent sum of 77(1440) and the broad CT; 
(b) shows the total J p = CT contribution; the lower histograms shows the broad CT intensity 
from Eq. (7) and the 77(1890) contribution; (c) the 1 + intensity; (d) experimental background 
(BG) and the 2~ contribution. 



intensity, also the broad component from Eqs. (7)-(9) and thirdly 77(1890). Panel (c) shows 
the small peak due to /i(1285). Panel (d) shows background as the upper histogram and the 
small 2" contribution required above 1600 MeV from t?(1645) and t?(1870). A detail is that it 
is necessary to check that the /o(980)t7 contribution to rj-KTi is consistent with data on rjK + K~ . 
The dashed curve on Fig. 1(f) shows the signal predicted for / (980)t7 with / (980) — > KK. It 
is consistent with data near 1600 MeV. 

The sharp rise of the r]irir mass distribution at 1495 MeV is not fitted perfectly in Fig. 
7(a): the data would prefer a rise 10 MeV lower than the fit. The rapid rise requires a strong 
contribution from / (980)t7. The fit which is shown uses the BES parametrisation of /o(980) [26]. 
A reparametrisation is desirable including the threshold cusp in the real part of the amplitude 
at the KK threshold. This has been tried, but does not cure the problem: it simply sharpens 
the / (980) peak near the KK threshold. 

Fig. 8 shows the fit to ^(KKn) data. The peak sits on a small experimental background 
shown in (a) and a slowly rising physics background, clearly visible above experimental back- 
ground. This physics background has J p = 0~ and is parametrised by the broad 0~ contribution 
decaying to K*K, nK and aon, do — * KK. The coupling constants of the first two are fitted 
freely; the last component is derived from the ao(980) fitted to rj-nn data using the branching 
ratio between rjir and KK in Ref. [24]. 

The KKtt peak of Fig. 8 is decidedly asymmetric. The main source of the asymmetry can 
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Figure 8: The fit to KKir data. In (a) the upper histogram shows the full fit and the lower one 
the J p = 1 + contribution from /i ( 1285), /i (1420) and /i(1510); the experimental background is 
labelled BG. (b) shows the 0~ contribution alone. 



immediately be traced to the s-dependence of the K*K decay channel whose contribution is 
shown in Fig. 8(b). The lower histogram in Fig. 8(a) shows the contribution from J p = 1 + . 
All fits require some 1 + contribution above /i(1420); this extra contribution is parametrised in 
the present fit by /i(1510) with mass and width quoted by the PDG. Changing the parameters 
to values of Aston et al. [27] changes log likelihood only by 3.6. 

Attempts have been made to fit the small shoulder in Fig. 8 at 1530 MeV in several ways. 
This mass is close to the opening of the KKrj channel at 1539 MeV. It raises the possibility of a 
drop in the line-shape associated with the opening of this channel; this is a unitarity effect - the 
opening of the r\KK channel robs KKtt of some signal. However, this is immediately excluded 
by the small branching ratio of r\KK compared to 7771-71". A more subtle possibility is that it 
may be associated with the /o(980)?7 channel, which also peaks at 1539 MeV, since /o(980) itself 
peaks at the KK threshold. The fit to 777771 data requires a strong contribution from /o(980)ry. 
However, it is spread out quite significantly according to phase space for the / (980)?7 final state, 
as shown below in Fig. 13(a); that is too broad to explain a narrow shoulder at 1530 MeV. Note 
that there is independent evidence for this / (980)7? contribution to J p = 0~ [28]. 

The / 2 (1565) — > [a 2 (1320)7r] L=1 does not account for the 1550 MeV structure in r/nn or the 
1530 MeV shoulder in KK7r; there is no visible a2(1320) peak at this mass in i]7T or KK and the 
angular dependence is wrong. Furthermore, the magnitude which can be fitted to ^(1565) in 
747T rules out the possibility of explaining the 1530 MeV shoulder with decays to [a 2 (1320)7r] /_,=!. 

One further possible explanation remains. In the accompanying paper on J/ty — > ^An, 
there is evidence that the pp peak at ~ 1600 MeV is actually resonant at ~ 1560 MeV. If this 
possible resonance has a weak decay to K*K, the shoulder at 1530 MeV can be fitted as an 
interference effect. The resulting fit is shown in Fig. 9 in identical format to Fig. 8; the addition 
of 7/(1560) with a width of 280 MeV improves log likelihood by 62. This explanation is clearly 
speculative. It is to be expected that two neighbouring resonances at 1440 and 1560 MeV will 
provide considerably more flexibility in the fit than 7/(1440) alone. However, the magnitude of 
the contribution from 7/(1560) is fairly small, as would be expected if its mixing with ss is small; 
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Figure 9: As Fig. 8 after the inclusion of 77(1565) . 



its intensity is 19% of that from 77 (1440). This additional resonance does not provide the 10 
MeV mass shift near 1500 MeV needed for a perfect fit to the r/Trn mass projection of Fig. 7(a). 
A possible remedy for this detail will be discussed in Section 7. 
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Figure 10: The fit to KKi\ data without /1 (1510). In (a) the upper histogram shows the full 
fit and the lower one the J p = 1 + contribution from /! (1285) and / 1 (1420); (b) shows the 0~ 
contribution alone. 



Let us now turn to the J p = 1 + contribution. Fig. 10(a) shows the poor fit without /i(1510). 
The remaining component due to /i(1285) and /i(1420) is shown as the lower histogram. The 
/i(1420) is parametrised including the s-dependence of the K*K channel and the associated 
dispersive term. The result is to extend the tail of /i(1420) to slightly higher masses, but 
does not substitute for /i(1510). The mass of / 1 (1420) is held fixed at the PDG value and its 
full width at half maximum is fixed to the width quoted by the PDG. The magnitude fitted 
to / 1 (1420) increases in Fig. 10 to 743 events compared with 488 in Fig. 8. The reason for 
this is the requirement for some significant 1 + signal in the mass range around 1500 MeV. 
An alternative is to increase the width of /i(1420) quite significantly. This possibility will be 
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discussed in Section 7. 



5.1 The fit with separate 7/(1405) and 7/(1475) 

The full s-dependence is fitted to both states including the dispersive terms. This requires a 
determination of the channels which contribute to each. 

The a (980)7r and r\a channels are naturally attributed to 7/(1405), as in earlier work. At- 
tempts to include these channels into 7/(1475) lead to negligible contributions. The K*K channel 
is naturally associated with 7/(1475), as in earlier work. If it is added to 7/(1405), log likelihood 
improves only by 3.9, which is insignificant. Here, it is necessary to comment on earlier work 
of Adams et al. [5]. Their Fig. 4(d) shows a low mass K*K component in addition to that of 
7/(1475). However, from the present analysis it is clear that this is an artefact arising from their 
free fit to a (980)7r, a KK. The resulting large a (980) intensity in KKtt is inconsistent with 
the r/n-K data by a factor 10, so a constraint on the magnitude of a (980) — > KK is essential. 
This eliminates any K*K contribution from 7/(1405). 

If kK is fitted freely to both 7/(1405) and 7/(1475), it goes almost entirely into 7/(1405). 
Removing the small component from 7/(1475) makes log likelihood worse by only 6.2, which is 
barely significant. Furthermore, it removes a destructive interference with K*K, symptomatic 
of instability. The final fit includes kK only in 7/(1405). 

The remaining component is /o(980)r/. It improves log likelihood of the fit by 162, a decisive 
amount. It plays a strong role in creating the sharp rise in the rjim mass projection of Fig. 
7(a) at 1500 MeV. One might therefore naturally associate it with 7/(1475). However, the fit 
contradicts this possibility. A free fit to both 7/(1405) and 7/(1475) leads to negligible contribution 
to 7/(1475) from / (980)t/. It is easy to locate the reason. Associated with the / (980)ry channel 
is a large dispersive term m(s) in the denominator of any resonance to which is couples. The 
phase variation created by this dispersive term is not compatible with 7/(1475) — > /o(980)t/. The 
final outcome is that 7/(1475) has only a K*K contribution and all other channels contribute 
via 7/(1405). 

The peak in r/Trn data is at 1415 ± 5 MeV. However, because of the interference with the 
broad 0~ contribution, the mass which can be fitted is quite flexible, with an uncertainty of ±20 
MeV. The fit actually optimises at 1429 MeV because the phase variation associated with m(s) 
in the Breit-Wigner denominator favours the higher mass. 

For 7/(1475), there is a strong correlation between the fitted width and the optimum mass. 
The strict optimum is at 1505 MeV. However, log likelihood is worse by only 4 if the PDG mass 
of 1476 MeV is used. In order to allow comparison with published results, the mass is fixed at the 
PDG average, 1476 MeV and the width is then optimised. The fitted width at half-maximum is 
78 MeV, compared with the PDG average of 85 ± 9 MeV. The fit to the r/mr mass distribution 
using these values is shown in Fig. 11. Near the 1415 MeV peak, it is somewhat different in 
shape to that of Fig. 7. The fit to the rapid rise at 1500 MeV is no better than in Fig. 7(a). 

The fit to KKn data is shown in Fig. 12. The mass distribution of the peak is fitted rather 
better than in Fig. 8. The reason is obvious: two resonances allow much greater flexibility. In 
particular, the phase variation with s allowed by 2 resonances has much greater freedom than 
the fixed relation between magnitude and phase for 7/(1440) alone. With two resonances, it is 
possible to arrange that the mass distribution in Fig. 12 rises rapidly and falls rapidly, but with 
something approaching a flat top. 
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Figure 11: The fit to the r/inr mass distribution using both 77(1405) and 77(1475). 



However, there is very little difference in log likelihood between fits with one resonance or 
two. The fit to 777T7T data with two states is better by 6.5 in log likelihood and the fit to KKn 
data is better by only 2.6, making a total of 9.1 for the two sets of data. This hardly warrants 
confidence in the need for two resonances. 




Table 1 shows changes in log likelihood for both types of fit when individual components 
are removed and magnitudes and phases of all other components are re-optimised. The third 
column shows results for fit A with one resonance, 77(1440), and the fourth shows results for fit B 
with two resonances. This Table also includes branching fractions for channels or combinations 
of channels (where interferences are large). A warning is that for r/(1440) there are strong 
interferences between kK and K*K and these are somewhat different for fits A and B. There 
are likewise large interferences between kK and K*K for the broad 0~ signal. Errors given for 
branching fractions cover uncertainties in interferences. Because of large interferences between 
the broad CT and 7/(1890), branching fractions of Table 1 do not add up to the total branching 
fraction given in Eq. (5). 

The one feature which might favour a fit with two resonances rather than one is the /i(1510) 
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contribution. Its parameters are fixed at PDG averages: M = 1518 MeV, T = 73 MeV. It 
improves the fit with a single r/(1440) and an /i(1510), which is regarded by some authors as 
questionable, by 54.0. The improvement for the fit with 2 resonances is 38.9. So the fit without 
any /i(1510) is better for fit B. Both fits require some /i(1510) contribution, or something 
resembling it. A test has been made of the requirement for /i(1510), as distinct from a CT 
contribution. A further 0~ state is added to fit B using r/(1405) and 77(1475) . It is given 
identical mass and width to /i(1510). The result is that the /i(1510) contribution dominates 
over that from 7/(1510), confirming a definite requirement for a 1 + component in this mass range. 
In Section 7, a possible explanation in terms of the line-shape of A (1420) will be presented. 
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Table 1: Changes in log likelihood when individual components are removed and all other 
components are re-optimised in magnitude and phase. Column 3 shows results fitting 7/(1440) 
alone and column 4 results fitting both 7/(1405) and 7/(1475). Columns 5 and 6 show branching 
fractions for J/\l/ — > 7 + X, multiplied by 10 4 . 
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Figure 13: (a) Contributions to the line-shape of ?7(1440), Fit A; (b) the variation with mass of 
phase space for each channel. The aovr signal refers to 777171; that for ao — > KK is 24% of the 
signal in r/nn. 

5.2 Further details of the fit with 7/(1440) 

Fig. 13(a) shows as the full curve the overall line-shape of ?7(1440), scaled down by a factor 
2 in order to display other contributions. The peak for K*K is at 1440 MeV, but there is a 
pronounced tail at high mass. One should note that in Obelix and Crystal Barrel data, this 
high mass tail is suppressed strongly by the phase space for production in pp annihilation at 
rest. 

From Fig. 13(a) it is quite difficult to decide on the centre of gravity of the K*K contribution. 
Averaging this contribution up to 1600 MeV it is at 1474 MeV, remarkably close to the PDG 
mass for 77 (1475), namely 1476 MeV. Putting this in reverse, if one unfolds the expected s- 
dependence of K*K, the resonance mass should be expected at 1440 MeV, and agrees with the 
fit with a single resonance. The /o(980) channel has a weak maximum near 1530 MeV, but also 
extends strongly to high mass. At high mass there are strong interferences between K*K and 
kK. Constructive interferences between them shifts the peak in KKir up to 1450 MeV. Fig. 
13(b) shows the dependence of phase space for each channel on mass, but with an arbitrary 
scale for each channel, so as to fit all the curves on to one figure. 

Fig. 14(a) shows the Argand diagram for the coherent sum of ao(980)7r amplitudes from 
77(1440) and the broad 0~ signal on which it sits. The diagram is drawn for an a (980) mass of 
980 MeV, so as to avoid the effect of phase space on the integrated amplitude. Numbers on the 
figure indicate rjirir masses in GeV. Fig. 14(b) shows the Argand diagram for the 77a amplitude 
evaluated for a a mass of 800 MeV. Up to 1415 MeV, the real parts of the amplitudes are large. 
The dispersive contribution to the real part is large because the amplitude senses the steep rise 
at higher masses in the imaginary part - an impending barrier, which gets closer as the mass 
rises. As the K*K channel opens, there are two effects. Firstly, it increases the magnitude of 
the Breit-Wigner denominator and damps all r/nn amplitudes. Secondly, the real part of the 
amplitude is approximately proportional to the gradient of the imaginary part. It turns over 
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Figure 14: Argand diagrams for (a) a (980)7r and (b) rja for the coherent sum of ?7(1440) and 
the broad CT contribution. Masses are marked in GeV. 



rapidly at 1420 MeV and plunges through zero on resonance at 1440 MeV; these two features 
produce the narrow peak at 1415 MeV in r/n-R. The / (980)7r amplitude contributes strongly 
to the steep rise in r/Tnr at 1500 MeV. The cusp at 1600 MeV in the broad 0~ contribution is 
largely responsible for the shoulder in r/TTn data at 1550 MeV. 

Table 2 gives branching ratios integrated from threshold up to 1.8 GeV. Systematic errors 
from this cut-off are a factor (1.0 ± 0.1) if the cut-off is reduced to 1650 MeV or increased to 
2100 MeV. Errors cover variations observed in many alternative fits; they also cover the (small) 
discrepancies between the fits with one resonance or two. Note that there is large interference 
between a (980)7r and rja. 



Channel 



Branching fraction (%) 



K*K 
kK 
rja 



40 ±5 
35 ±5 
4.1 ± 1.6 
5.5 ±2.1 
1.3 ±0.5 
14 ±3 



a n,a 
a 7r,a ■ 



KK 



W80)rj 



Table 2: Branching fractions of ^(1440). 



5.3 The branching ratio KKh/tjittt 



From Table 2, the ratio 

BR(rj(UA0) -> KKtt) 
BR(r](lU0) -> 777T7T) 




(12) 
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Figure 15: The factor H{m) allowing for the phase space and dynamics of the production process 
pp r7(1440)<7 as a function of the mass m of 7/(1440). The full curve shows production via 
the a pole for the spectator pair and the dashed curve the result using the tttt elastic scattering 
amplitude. It is normalised to 1 at 1440 MeV. 



The errors have been obtained by adding in quadrature the errors for K*K and ktt to get overall 
branching fractions for KKn in the range (82 — 68)%, then dividing by the corresponding r]Tm 
fraction. 

In their review, Masoni et al. [8] present in their Table 8 their estimate of the same branching 
ratio, or more exactly its inverse. Converting their entries to the definition used here, they 
estimate values in the range 0.98 ± 0.37 from Mark III and DM2 data, 0.30 ±0.1 using Crystal 
Barrel and bubble chamber data of Baillon et al. [29], or 0.11 ± 0.05 from Crystal Barrel and 
Obelix data. They use the disagreement between these sources as an argument in favour of 
separate 77(1405) and 77(1475) . 

However, this analysis is misleading. When one analyses data for pp annihilation, it is vital 
to allow for the suppression of the K*K channel by the available phase space. This was ignored 
by Masoni et al. The factor H(m) multiplying the line-shape of r/(1440) is: 




Here Si is the mass squared of the spectator txtx pair, and p pro dn is the momentum of ?7(1440) or 
the spectator a in the pp rest frame. Next, M p and m n are the masses of the proton and pion. 
Also T a is the scattering amplitude for the spectator txtx S-wave pair and F is the form factor 
for the production process, taken here as exp(— 2.25p pro d n ) with p pro d n m GeV/c. 

The variation of H{m) with the mass of 7/(1440) is shown in Fig. 15. It is not clear whether 
the nn spectator pair will be produced via the a pole or via the nn elastic scattering amplitude. 
The former is shown as the full curve and the latter as the dashed curve. They are sufficiently 
close that one can work with the average. The high mass tail of 7/(1440) is cut off strongly in 
pp annihilation. It goes to zero at 1.60 GeV: the highest mass which can be produced in pp 
annihilation with two spectator pions. 

There are two consequences. Firstly, the factor H{m) amplifies the lower side of 7/(1440) 
with respect to the upper side. For the narrow 77a" and a (980) channels, the effect is to shift the 
apparent peak of 7/(1440) downwards by 4.7 MeV; i.e. the Breit-Wigner mass is actually 4.7 MeV 
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higher than is fitted with a Breit-Wigner amplitude of constant width. For the kK channel, 
the effect is larger because the phase space for this channel, shown on Fig. 13(b), increases 
significantly with m. The result is that the mass fitted to this channel is moved downwards 
by 9.7 MeV, so the true mass fitted as a Breit-Wigner resonance of constant width should be 
increased by this amount. 

The second effect is that a considerable part of the KKtt signal is lost in data on pp anni- 
hilation; all of it is lost above 1.60 GeV. 

A check will be presented here on the ratio KKn/rjim from data on pp annihilation. Atten- 
tion will be drawn to many sources of uncertainty. The check is done using data of Anisovich 
et al. [6] for pp — > (i]tt + tt~)(tt + tt~) and bubble chamber data of Baillon et al [29]. Table 2 of 
Anisovich et al. gives branching fractions for all amplitudes fitted to their data, in particular 
77(1440) — > rja and ao(980)7r. This Table reveals an important point. There is large destructive 
interference in those data between these two amplitudes. This destructive interference has a 
strong effect on the branching fraction quoted for 7/(1440). However, in reality there are major 
uncertainties in determining the interference. There are actually four charge cominations in 
those data allowing considerable freedom in interferences. The individual contribution from rja 
is quoted as 2.31 x 10~ 3 , that from a (980)7r is 9.6 x 10~ 3 and the interference between them 
contributes —2.28 x 10~ 3 . It clearly makes a large difference how this interference is taken into 
account. 

Relative magnitudes of all decay channels of a resonance should be the same in all production 
processes. The procedure adopted by Anisovich et al., and used here, is to derive coupling 
constants from the fit to data, then evaluate the intensity of amplitudes by taking the modulus 
squared of individual amplitudes. Interference between decay channels of one resonance are 
kept, since these are a property of the resonance. Interferences with other resonances in the 
data are ignored because such interferences are a feature of the particular environment in which 
the 77(1440) is produced. 

The resulting branching fraction for production of 7/(1440) in i]tttt via pp annihilation is 
1.09 x 10~ 3 in the observed charge states. This needs to be multiplied by a factor 9/4 to 
allow for other charge states, giving a result of 2.45 x 10~ 3 for pp — > 7/(1440), 7/(1440) — > t]tttt. 
The result for all charges of KKtt quoted by Baillon et al. is (2.0 ± 0.2) x 10~ 3 . This gives 
a ratio #1440 = BR(r}(lU0) -> KKn)/BR(r}(U40) -> t/tttt = 0.82. Masoni et al. use the 
smaller branching ratio 0.61 ± 0.19 found by Amsler et al. [30] in an earlier analysis of pp data. 
However, that publication does not specify if any allowance is made for the effect of destructive 
interference between 7/(1440) — > r/a and a (980)7r. 

The effect of H{m) on the integrated intensities for final states rja and ao(980)7r is small 
because the phase space p(s) for these channels varies slowly with m and there is an approximate 
cancellation between upper and lower sides of the resonance. For the kK channel, there is a 
noticeable effect because of its fairly rapidly increasing phase space, see Fig. 13(b). Using Fig. 
15, the kK component in pp annihilation needs to be multiplied by a factor 1.38 and the K*K 
contribution needs to be multiplied by ~ 3.1. If we take the relative contributions of K*K and 
kK to be in the ratio 40:35 of Table 2, the result is to increase the KKtt branching fraction of 
Baillon et al by a factor 2.3, resulting in a value -R1440 of 1.88. 

This is still a little smaller than the value derived from Table 2, namely 3.0^o!g- However, the 
essential conclusion is that the branching fraction between KKtt and t\tttt is subject to a large 
systematic correction when evaluated from pp data. The result from J/\I/ data does not suffer 
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from these systematic problems. Note that the line-shapes fitted to present data are perturbed 
by the slowly varying factor for production with L' — 1; this factor and the associated form 
factor have been included in the fit and the arithmetic of branching ratios. One should also 
note that the signal assigned by Obelix to the KKir channel is not visible as a peak in the data, 
but is obtained from the difference in the line-shape of the KKir signal from a Breit-Wigner 
of constant width fitted to 7/(1405). A further detail is that the phases of individual decay 
channels are affected by multiple scattering with other components of the fit. They may differ 
from one production process to another. As a result, it is difficult to give a definitive result for 
the branching ratio between KKn and r/nn. 

The reality of the situation is that there is major flexibility in the magnitudes of contribution 
which can be fitted to 7/(1440) in rjim channels. In the data presented here, there is flexibility in 
what can be fitted to the broad interfering component. This flexibility is apparent from the fact 
that fits with a single 7/(1440) or separate 7/(1410) and 7/(1475) differ very little in log likelihood. 
There is yet further flexibility if the r\a channel receives additional contributions proportional 
to the nn elastic amplitude. 

6 Fits to data above 1600 MeV 

In i]7ttt data shown in Fig. 7, the broad 0~ signal in Fig. 7(b) obviously fails to explain fully 
the 7]7ttt signal above 1600 MeV. The K*K* threshold may contribute something to the drop 
above 1870 MeV, but alone does not fit the high mass spectrum. The remedy adopted here is to 
fit a resonance form as an expedient, although there is no phase information, hence no definite 
evidence for a resonance. The optimum fitted mass and width are 1890 ± 20 MeV, V = 260 ± 35 
MeV. The strongest decays are to ao(980)7r, but there is also a rather strong decay to /o(980)t/. 
Fig. 16 show mass projections from r]7i7i Dalitz plots on to ixix and r/n mass. The / (980) is 
clearly visible in Fig. 16(g) at 1 GeV. In Figs. 16(h), there is a a small but definite a 2 (1320) 
peak. It does not fit naturally to decays of 7/ 2 (1870) and/or 7/ 2 (1645). It is of roughly the right 
magnitude to originate from a 2 (1320)p background. The p(770) signal from this background is 
also just visible in Fig. 16(g). 

On the upper side of the 1850 MeV peak is a small 2~ component, which may be fitted 
with 7/2(1870). When this contribution is omitted from the fit, log likelihood gets worse by 84. 
The relative intensities in a (980)7r, a 2 (1320)7r, / 2 (1270)t/ and arj decays are fixed to accurate 
values from Crystal Barrel data [31]. Production phases, attached to each amplitude in the 
isobar model, are fitted freely. So there are five fitted parameters. Statistically, the significance 
of the 2~ component is 10 standard deviations, despite the signal being small. The reason 
for this significance is that amplitudes have a distinctive angular dependence which is readily 
identified. For completeness, a contribution is also included for 7/ 2 (1645). It is much smaller and 
barely significant. In Ref. [31], it is suggested that 7/ 2 (1645) is a standard qq state, the isospin 
partner to 7r 2 (1670); t/ 2 (1870) is interpreted as a hybrid. This interpretation is consistent with 
the observations here. The production of a hybrid is naturally explained in terms of cc — > ^ygg, 
followed by conversion of one gluon g to qq. 
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Figure 16: Fits to mass projections on to (a) M(itit) and (b) M(nr]) for the rjiriT mass range 
1230-1330 MeV, compared with the fit (histogram); (c)-(d) the mass range 1350-1480 MeV; 
(e)-(f) 1480-1600 MeV; (g)-(h) 1700-1900 MeV. 
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6.1 The magnitude possible for X(1835) — > r]7T7T 

The BES Collaboration has presented evidence for an X(1835) in n'inr with a production branch- 
ing fraction BR(J/m -> 1 X).BR(X -> tt+tt"^') = [2.2 ± 0.4(stat) ± 0.4(sys£)] x 10" 4 and a 
width of 68 MeV [11]. There is no obvious signal in 7/7r + 7r~ data of Fig. 7, although the fit to 
data is slightly ragged in this mass range. 

A free fit including it with published mass and width gives 189 / (980)7T events, 94 a (980)7r, 
164 no and 86 a 2 (1320)7r. Log likelihood improves by 58 for 8 additional fitting parameters. 
There is flexibility in the fit, and the total number of events can be increased from 533 to 1182 
if log likelihood is allowed to get worse by 2 for every channel, corresponding to an overall 
4 standard deviation change. This is a cautious limit taking account of the fact that the 
components fitting this mass range are uncertain. The corresponding branching fraction is 

n ) = 0.38 x 1(T 4 . No statement can be made about its 
possible coupling to KKn, since the mass range of the present data extends only to 1800 MeV. 

Huang and Zhu [32] point out that its mass, total width, production rate and decay pattern 
favour its assignment as the second radial excitation of n'. More precisely, the small branching 
fraction in nnn favours a mixing angle between nn and ss close to that of the rj'. The rj' (958), 
77(1440) and X(1835) would then form a Regge trajectory with a spacing in M 2 of ~ 1.22 GeV 2 , 
close to that of other well establshed trajectories [33]. The peak observed by BES2 at 2240 MeV 
in 00 would fit naturally as the fourth member of this sequence [34]. However, resonant phase 
variation needs to be established for 7/(1835) and the 00 peak. 

Ma [35] tries to interpret 5(1835) as a baryonium state associated with the threshold peak 
observed in pp in BES data [36]. This interpretation leads to a decay stronger to nun than n'nn, 
and is clearly ruled out by the present data. 

7 Interpretation of 7/(1440) and /i(1510) 

If there is just one 7/(1440) in the mass range 1400-1500 MeV, it either has a large SU(3) 
singlet qq component or a glueball component. Either would explain its strong production in 
J/ty radiative decays, where cc appear in a singlet combination. It could likewise be produced 
strongly in pp annihilation. However the K*K channel is SU3 octet, so the 77(1440) (or 77 (1475)) 
cannot be purely singlet. 

It is well known that isolated ss states such as / 2 (1525) are produced weakly in pp annihila- 
tion. If the 7/(1475) is to be produced strongly in pp annihilation, it should have a substantial 
nn component, but that is inconsistent with decays purely to K*K. The 7/(1405) decays much 
more strongly to kK and / (980)ry than to no + a (980)7r. This would not arise naturally if it 
is a pure nn state. So the conclusion is that a substantial ss (or gluball) component must be 
present in 7/(1440) or 7/(1405) +7/(1475). 

The existence of 7/(1295) depends primarily on data in three papers on n~p — > (nim)n and 
one on ir~p — > (KKn)n [1]. In this process, there should be a clear distinction between 77(1295) 
which decays via the S-wave and /i(1285) which decays via the P-wave. If production is limited 
to a single exchange process, these are immediately distinguishable. In pp annihilation the 
situation is much more complex, since there are four combinations of nmr produced in every 
event [6], making the separation between 7/(1295) and /i(1285) much more difficult. Those 
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reviews which argue against the existence of 77(1295) base their arguments on the fact that 
there could be more than one exchange mechanism in 7r~p production data, leading to possible 
confusion between the flat angular distribution from 77(1295) and rank 2 processes involving 
/i(1285). It is unfortunate that the papers of primary importance published on r/(1295) do not 
give quantitative estimates of its significance. If the data still survive, it would be valuable if 
the groups concerned would publish values of log likelihood or x 2 showing quantitatively how 
their fit changes when the 77(1295) is omitted. This information would sharpen the arguments 
for or against 77(1295) . 

In the data presented here, most of the peak near 1285 MeV certainly is due to /i(1285); 
a free fit including 77(1295) improves log likelihood only by 3.7 for four parameters fitted to 
magnitudes and phases of decay amplitudes to ao(980)7r and rja. Seven events are fitted to 
77(1295) compared with 230 to /i(1285). So the evidence for r/(1295) in present data is not 
significant. 

Is the 7/(1440) really split into two components? Although fit A optimises 77(1405) at 1429 
MeV, a persistent feature of data fitted to 77(1405) in other analyses is a mass 1405 MeV in i]7iit 
and 1414 MeV in kK . As shown in Section 5.3, data from pp annihilation require corrections to 
these masses for the phase space in the production process. This shifts the mass fitted to r\ixix 
data up by 4.7 MeV and that in KKtc data by 9.7 MeV. Weighting these with the intensities 
fitted to r/TT-K and kK in Table 2, the mean mass of 7/(1405) becomes 1421 MeV. There is then 
a discrepancy in mass of 18 ± 5 MeV with the mass in fit A with a single resonance. This is not 
convincing evidence for two separate resonances near 1440 MeV. However, it must be said that 
one cannot move the 7/(1440) much below 1439 MeV without a poor fit to the peak in KKn. 
How could the difference of 18 MeV be explained between rjirn data and KK? 

Klempt [37] has pointed out a possible explanation. A radial excitation has a node in its 
wave function causing a zero in decays to a n in this general mass range. Firstly this attenuates 
the 7/7T7r signal and may explain why it is weaker than kK and K*K. Secondly, such a node 
introduces a further s-dependence into rjim. Suppose the amplitude for the a 7r channel is 
multiplied by a factor (s — s)/s , where s is the location of the zero in the amplitude. A 
resonance at 1439 MeV can have a peak at 1410 MeV if there is a zero in the amplitude at 1.47 
GeV, just where there is a dip in the r/im mass spectrum, Fig. 7. However, the peak is attenuated 
more on its upper side than its lower side and the peak is somewhat asymmetric. Tests on the 
present data show that this explanation does cure the 10 MeV discrepancy between tj-kti data and 
the fit near 1500 MeV. However, a factor 10 higher statistics are needed to test the idea precisely. 
It cannot be tested in pp annihilation because of the four different charge combinations. Zeros 
would appear at different masses in r\a and kK, depending on matrix elements, leading to a 
potentially complicated line-shape with peaks at different masses in different channels. 

There is one further possible complication. The dominant decay of 7/(1440) is through KKn. 
There is inevitably some rescattering of KK through ao(980) to the i]7itt channel via what are 
called triangle graphs. This was discussed in some detail in Ref. [6]. It will have the effect of 
increasing the KKn decay branching fraction of 7/(1440). 

Let us turn to 77 (1890). The general pattern of the broad components in the 7/7T7t mass 
projection of Fig. 7 is that there are two components peaking at 1890 and 1600 MeV. It is 
possible that the latter is resonant. In the jAtt data there is evidence for a resonance at 1970 
MeV with a width of 210 MeV. If the parameters of 7/(1890) are changed to those values, log 
likelihood for rjim data is worse by 42. This is a modest amount, in view of the systematic 
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uncertainties about how to fit the r/TriT mass range above 1600 MeV and the strong K*K* 
threshold nearby. It is therefore possible that these two signals arise from the same resonance, 
though this cannot be proved at the moment. 

Consider next /i(1420) and /i(1510). In Fig. 8(a), one sees that the effect of /i(1510) is to 
produce a slight bump at 1490 MeV and extend the tail of /i(1420). Longacre [38] has suggested 
that /i(1420) may be an abnormal state (or 'molecule' in current terminology) with an L = 1 
pion orbiting a KK S-wave combination. In this picture, there are final state interactions from 
two combinations of K* and a further final-state interaction between the two S-wave kaons. 
My suggestion is that these final state interactions enhance this state strongly at the lowest 
KK masses, but the rapid rise of the K*K phase space with momentum produces a long tail 
extending into the mass region of what is called /i(1510). 

It must be said in favour of the data of Aston et al. that they identified the foi(1370) 
J PC = 1~ + state and it was subsequently confirmed by Crystal Barrel data [39]. Furthermore, 
there is a 2~ + peak in their data which they did not claim as a resonance, but was subsequently 
identified as a resonance by Crystal Barrel [40] and confirmed by WA102 [41]. The conclusion 
is that the data of Aston et al. appear reliable. The primary question is one of interpretation 
of their J PC = 1 ++ signal. They did not include /i(1420) in their fit. So it is possible that the 
/i(1420) has a more extended line-shape than has previously been fitted. 

8 Summary 

The main experimental results are the observation of three peaks in the rjini mass spectrum: a 
narrow peak at 1415 MeV, a shoulder at 1550 MeV and a broad peak at ~ 1850 MeV. In KK7r, 
there is an asymmetric peak at ~ 1450 MeV. The asymmetry is due to the s-dependence of the 
K*K channel; its Dalitz plot immediately identifies a dominant K*K decay. The intensity for 
that decay increases rapidly from its threshold at ~ 1390 MeV as A; 3 , where k is the momentum 
of decay K and K* in the K*K rest frame. 

An essential point from the present analysis is that one must include this s-dependence 
correctly into the Breit-Wigner amplitude. This reduces the mass of 77(1475) to 1440 MeV. 

When analysing pp data, it is necessary to allow for the variation of phase space for produc- 
tion of rinn and kK with mass. This corrects the mass of 7/(1405) upwards by 4.7 MeV in r\ixix 
and by 9.7 MeV in kK. As a result, the mean mass difference between the possible 77(1405) and 
771440) is no larger than 18 ± 5 MeV. 

There is a dominant broad 0~ amplitude in rjiriT. When this is included into the partial wave 
fit, interference with a single 7/(1440) can account naturally for the observed peaks at 1415 MeV 
in r/TT-K and at 1450 MeV in KKir. There is no evidence for two separate 7/(1405) and 7/(1475) 
from the present data. 

The 7/(1835) is not observed in the present 7]7m data, and an upper limit has been derived 
for the branching fraction for its production in J/ty decays. This limit is a factor 5.7 ± l(stat) ± 
l(syst) smaller than the 7/(1835) observed in the 7/7T7r channel. The natural explanation is that 
7/(1835) is a radial excitation of 7/' and 7/(1440). 

As an aid to further work, the fully annotated Fortran code for the broad 0~ component and 
7/(1440) is available from the author. This code includes numerical values of g 2 from the present 
work, but they would need to be re-optimised if further data become available. 
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